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Abstract 

We determine the charm-quark mass mc{mc) in the MS scheme using measurements of charm 
production in deep-inelastic ep scattering at HERA in the kinematic range of photon virtuality 
5 GeV^ <Q^ < 1000 GeV^ and Bjorken scaling variable 10^ <x < 10 ^ The extraction of 
mc from this process with space-like kinematics provides complementary information to results 
from hadronic processes. The QCD analysis of the HERA data yields a value of mc{mc) = 1.27 it 
0.05 (exp)+j!;j!J (scale) GeV at next-to-leading order and of mc(mc) = 1.36±0.04(exp)l!5;j!J(scale) ± 
0. 1 (theory) GeV at approximate next-to-next-to-leading order. The results are consistent with and 
of comparable precision as the world average. 



'Permanent address: DESY, Notkestrafie 85, D-22607 Hamburg, Germany 



1 Introduction 



Quark masses are fundamental parameters of the Standard Model. However, due to confinement 
no free quarks are observed in nature. Therefore, for the determination of heavy quark masses 
a careful theoretical description is needed for appropriate observables. This then enables the de- 
termination of heavy quark masses by comparing quark mass dependent theoretical predictions 
to experimental data. In such an analyses reference must be made to the specific scheme used to 
define the quark mass and it is mandatory to include radiative corrections in Quantum Chromody- 
namics (QCD) beyond the leading order (LO). 

In QCD predictions for hard scattering processes employing the pole (or on-shell) scheme for 
the definition of the heavy-quark mass, ntq is chosen to coincide with the pole of the heavy-quark 
propagator at each order in perturbative QCD. It is known since long, that the concept of pole 
masses in QCD has an inherent drawback, see, e.g., Q. Since quarks are confined inside hadrons, 
there is no pole in the quark propagator of the full theory. As a consequence, the quark mass 
parameter in the pole mass scheme is limited to the perturbative domain with corrections of order 
0{A.QCD/mq) and its value depends strongly on the order of perturbation theory. Alternatively, 
QCD predictions for the charm production cross section in DIS can be considered as a function 
of the running mass mg(/j, ) at the renormalisation scale /j^ in the MS scheme. The quark mass, 
nigdj,-), is treated on the same footing as the strong coupling constant (Xs{iUr) and represents a 
prominent example of the so-called short-distance mass definition, which requires the heavy-quark 
masses to be evaluated at the scale /Ur much larger than the QCD scale Aqcd, i-e-, /Jr ^ Aqcd- 

To this date, the charm-quark mass has been measured very accurately from data in electron- 
positron annihilation, e.g., with the help of QCD sum rules, or from numerical simulations on 
the lattice, see HZj. These determinations are either entirely non-perturbative (lattice) or related to 
scattering processes in time-like kinematics only. It is therefore of particular interest to examine the 
possibility of alternative charm-quark mass determination from heavy-quark production in deep- 
inelastic scattering (DIS). This reaction proceeds in the space-like kinematics with the distinct 
structure of the QCD amplitudes and in this way it provides an additional check of the QCD 
parameters' universality. 

In this paper the theoretical framework developed in [j3| is applied to determine the charm- 
quark mass at next-to-leading order (NLO) and approximate next-to-next-to-leading order (NNLO). 
The recent results from the HI experiment (U-O on charm production in deep-inelastic ep scat- 
tering at HERA are used in addition to the data sets of [|71- By using the MS scheme the current 
analysis profits from improved stability of the perturbative expansion and the reduced theoretical 
uncertainty due to missing higher order contributions estimated by the variations of the renormali- 
sation and factorisation scales. 

This paper is organised as follows: after a brief reminder of the theoretical basis for the treat- 
ment of heavy quarks in DIS we present the determination of the charm-quark mass in the MS 
scheme, based on recent charm production measurements from HI. Special emphasis is put on the 
analysis procedure to account for the charm-quark mass dependence of the D*^ production cross 
section in the fiducial kinematic range. Finally, the running charm-quark mass md/Jr) is deter- 
mined at NLO and at NNLO in QCD. In the NNLO case, the theory predictions are available only 
to an approximation with a substantial uncertainty arising from missing information in the NNLO 
DIS Wilson coefficients. The results are compared to previous extractions of mc{mc) from the DIS 
data in |l3l[8]| and to the results of different determination methods, used for the world average [|2]|. 
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2 Theory framework 



Charm-quark production in neutral-current DIS proceeds by scattering of a charged lepton e off a 
proton P following the reaction 

e{l) -f P{p) ^ e(Z') c(pi) c-(;?2) + ^ , (D 

in which a virtual boson of space-like momentum = —q^ = — (Z — l')^ is exchanged. Here 
/, p, I', pi and p2 denote the four- vectors of the incoming electron and proton and the outgoing 
electron, charm quark and anti-charm quark, respectively. The charm-quark pair, cc, in the final 
state is heavy, such that in^. ^ Aq^q holds. We restrict the momentum transfer to be much smaller 
than the Z-boson mass 

The inclusive cross section is expressed in terms of the standard DIS proton structure functions 
Fk with k = 2,L. The charm contributions to the inclusive structure functions are denoted 
by F^^{x,Q^,m^). In referring to reaction ([T]) the structure function Fj^^ requires by definition a 
charm-quark pair in the final state, but it is otherwise a completely inclusive quantity, especially 
with respect to the proton initial state. Thus, the cross section for the reaction ^ reads 

^ = ^Kl+(l-3')'}^f>,e^-?)-/^r>,e^-?)] • (2) 

Here a is the electromagnetic coupling constant and the DIS variables x and y are defined as 
X = I {2p ■ q) and y = {p ■ q)/{p- 1), respectively. For later use the reduced cross section is 
introduced as 



GUx,Q\ml) = Ff{x,Q\ml) - -—^-—^F[' [x^Q^m]) . (3) 

In the standard factorisation approach to perturbative QCD the structure functions Fj^'^ can 
be written as a convolution of parton distribution functions (PDFs) and Wilson coefficients, see 
e.g., illlQl, 

Ff (-,e^m?) = ^ __E ^ /' f /, Q„(U,^},;.?) . (4) 

where z^^^ =1/(1 +4m^/2^) and e^. = 2/3 is the normalised charm-quark charge. The PDFs 
for the parton of flavor i are denoted as fi{x,ijj) and the sum in equation (|4]) runs over all flavor 
combinations, i.e., singlet and non-singlet, and the gluon. The Wilson coefficients Q., depend on 
the kinematic variables z and ^, 

Q-^ + s 

with s denoting the partonic centre-of-mass energy. 

For the treatment of heavy quarks in DIS the so-called fixed flavour number scheme (FFNS) is 
chosen with the number of active quarks in the proton rif = 3 and massive charm quarks appearing 
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exclusively in the final state. Moreover, the strong coupling constant is defined in the same scheme 
as (Xs{nf) with rif = 3 for charm quark production. This description of QCD with one massive and 
Hf light quarks can be related to QCD with (nf + l) light quarks by means of the standard matching 
conditions, cf. e.g., [[TT1IT2| . and for a discussion of all-order res ummations of logarithms in Q^/ni^ 
in so-called variable flavor number schemes (VFNS) llT3l . see [[141 [T5l . The Q , are obtained in 
perturbative QCD as an expansion in = CLs{fir), 

oo oo y 2 

CkAz,^y) = I {4iia,y4j{zA,P^) = £ {4iia,y £4f (z,^) In^^ , (6) 

here the renormalisation and factorisation scales are set to /j = fjf = /jy. In order to estimate the 
uncertainties due to missing higher orders, the scale /u is varied by a factor two up and down around 
the nominal value El 



/J e K^Q^ + Ainl with ^ < k < 2. (7) 

Conventionally, the Wilson coefficients are presented using the pole mass scheme for m^-, see 
e.g., [IHIlll. The necessary formulae for the conversion to the running mass mc(/Jr) in the MS 
scheme are well known [[T7l - [T9|[ and the application of the MS scheme to the calculation of ^1'*^' has 
been detailed in [|3|[ (see also [[20l[2TI for related work on heavy-quark hadro-production). For the 
inclusive cross sections at short distances an appropriate choice for the scale of the running mass 
mc{fJr) is = mc. The renormalisation group evolution for the scale dependence governed by the 
corresponding quark mass anomalous dimension [|22l[23l converges even for scales as low as the 
charm-quark mass, that is 0(1) GeV. 

The Wilson coefficients in equation ^ have been computed exactly to NLO and the functions 
cfV [0 are often used via the parameterisations of [[TOl , see also [|24l . At NNLO approximate 

(2\ (2') 

results for the most important gluon and quark coefficient functions, ^ and ^ , are known [[T6][ 
and denoted by NNLOapprox in the following. These are based on recent partial NNLO improve- 
ments for the structure function ¥2^^ which encompass various kinematic limits: 

(i) Cjg^*"^ near the partonic threshold for 5 ~ 4m^ , 

(ii) Cj^j ^™ and c^]^ at asymptotically high scales ^ , 

(iii) cPj™^^' ''andcPj™^^' at high energies 5 (small-x). 

The threshold approximation c'^^^ has been determined to the next-to-next-to-leading logarithm 
(NNLL) (see [|25ll26l for previous approximations at the level of the next-to-leading logarithm 
(NLL)). The function Cj^^*'^ is consistent with zero to the accuracy considered. Likewise, fully 

analytic results for / , i = g, q, in the asymptotic regime of Q ^ have been obtained. 
The corresponding calculations make use of the formalism of [[m[T2|[ and a number of lowest 
even-integer Mellin moments for the necessary heavy-quark operator matrix elements at three 



Most of current global PDF determinations assume /j, = /j/ = 2 in fits to DIS data. This is the appropriate scale 
choice for massless structure functions and in general for large values of Q when mass effects are negUgible. 
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loops ll27l - |29l . see also ll30l . In the high-energy limit the expression for ■ ™ ^ is exact to 
leading-logarithmic (LL) accuracy at small-x due to II31II but is approximate only at NLL. 

The combination of all available information leads to an expression for the NNLO Wilson 

(2) 

coefficients, c\ / of the form: 

4? - 4?*^+(l-/(^))cg^^■"+/(^)cg?^■"^"-^ (8) 

with a suitable matching function /(^) connecting the regions ^ and s ^ in^., cf. equa- 
tion (4.9) in lfT6]l . 

The approach chosen has one caveat, though, which is a non-negligible theoretical uncertainty 
due to the poorly constrained NLL term at small-x in c^^ ^ and the limited number of known 

three-loop Mellin moments [|27l - |29ll . To account for these deficits two different scenarios c^j'^ and 

(2) B ' 

c\ • in equation ([8]), carefully designed to parametrise the related uncertainties, will be considered 

in the subsequent analysis. For the exact definitions of c'^}''^ and c'^}'^ see equations (4.17), (4.18) 
and equations (4.21), (4.22) in lfT6ll . An additional theoretical uncertainty on F^*^ is estimated by 
the variation of the renormalisation and factorisation scales. This part of the theoretical uncertainty 
is obtained using exact results since the ^u-dependence is fully known at NNLO [|3l|25]| from the 
renormalisation group, see also lfT6l . 



3 QCD analysis of charm production measurements 

To study the impact of the charm production measurements at HERA on the determination of the 
MS charm mass a variant of the ABMl 1 fit is performed with the data of [l4]-[6l included. The 
data were collected by the HI experiment during the HERA II running period corresponding to 
an integrated luminosity of about 350 pb^^ and cover the kinematic range of photon virtualities 
5 GeV^ < < 2000 GeV^. The charm mass mc{mc) is a free parameter of the fit comparing 
equation ([3]) to the experimental data. Experimentally charm-quark production in DIS is tagged 
via fully reconstructed charmed mesons or by using secondary vertex information of tracks, ex- 
ploiting the longevity and the large mass of charmed hadrons. These measurements are restricted 
(e.g. in the transverse momenta or angles of the produced particles) by the acceptance of the de- 
tector. The phase space, in which charmed hadrons can be fully reconstructed is usually referred 
to as the visible or fiducial phase space. Corrections for the non-measurable phase space and the 
fragmentation of charm quarks to charmed hadrons are applied in the fit. In the following, the mea- 
surements of charm production used as input and their treatment in the QCD analysis are described 
with particular emphasis on the corrections for non-measurable phase space. 

The c-quark production in |l4l[5]| is tagged via fully reconstructed D*^-mesons in the decay 
mode D*^ (D^ ^^71^)71^. These measurements are restricted to the visible kinematic range 
of the D*^-meson's transverse momentum pj{D*) > 1.25 GeV and pseudo-rapidities |ri(Z)*)| < 
1.8 at medium virtuaUties < 100 GeV^. At high virtuaUties 100 GeV^ < < 1000 GeV^ 
the visible range of pt{D*) > 1.5 GeV and |ri(D*)| < 1.5 is covered. Due to the phase space 
limitations equation ^ cannot be used directly in the analysis of these data. Therefore, for these 
measurements the invisible phase space region is accounted for in equation Q through a factor 

£vis = Ovis(Z)*^)/a^^n, (9) 



4 



where (5^,is{D*^) and a^^^ are the QCD predictions for the D*^ meson cross section in the visible 
phase space and the charm production cross sections in the full phase space, respectively. These 
predictions are calculated at NLO in the FFNS with the fully exclusive program HVQDIS [[24||. 
The contribution from Z?-quarks to the inclusive D*"^ -meson production cross section, reported 
in EllSl, is estimated using HVQDIS and is subtracted. The Z)*^ cross sections are re-calculated 
using the recent branching ratio values [2]. 

The value of Eyis depends strongly on and its dependence on is taken into account in the 
fit. For this purpose Eyis is calculated in a first step for selected values of the charm quark mass, 
which correspond to the range of mc{mc) scanned in our fit: 



nicimc 



0.9, 1.05, 1.2, 1.35 GeV. 



(10) 



Since HVQDIS is based on the pole-mass definition, the calculations are performed using the 
following values for the charm quark mass in pole definition: 



m 



P"!'^ = 1.41,1.53,1.67,1.81 GeV, 



(11) 



obtained by mapping of the running-mass grid in equation (ITOI) according to the matching between 
pole and MS definitions at the appropriate order as detailed below. The relation is known to 

o(a3) m\m. 



m 



pole 



mAMr 



1 + -a + (-3. 1242 + 13.4434) + 1 16.504a3 



(12) 



where a = as{mc{mc)) /n. 

In the calculation of Eyis, the proton structure is described by special PDF sets, provided for 
this analysis. They correspond to the NLO variant of ABM 11 fit Q and are performed in the 
MS definition of charm quark mass. The PDFs depend on the value of through substantial 
correlations between a^, and the gluon PDF. In order to provide a fully consistent treatment 
of the charm-mass effects in Eyis one has to take into account this dependence. Therefore we use 
the PDFs, which exactly correspond to the current value of appearing in the fit. These PDFs 
are given by interpolation between those obtained in the variants of ABM fit with the m-c settings 
of equation (flOl) . The full dependence of Eyis on rtic is provided by parabolic interpolation of the 
HVQDIS results obtained at the values of rtic in equation (fTTI) with the interpolation coefficients P, 
calculated independently for the LO and NLO terms in CJ^-j,^^^jj, 



rc,LO I ^2 



pole -J 



^cc,NLO / y^2 
^vis/full (^'2 '3^'"^ 



pole\ 

c ) 



a?! 

i=0 



1,!- 



s/full 



(13) 
(14) 



Such a representation allows the determination of Eyis in terms of the MS mass. This is achieved by 
substituting n^f^^ in equations (fT3l) and (fT4l) with the matching condition of equation (fT2l) similarly 
to the approach used earlier to derive the heavy-quark electro-production coefficient functions in 
the MS mass definition The terms of O(a^) and higher appearing after substitution exceed the 
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NLO accuracy therefore they are dropped and the final expressions for CJ^^^^f^jj employed in our 
analysis read 



with 



vis /full/ ^2 / 

vis/full / „2 f \\ 



pVis/full 
M,0 

pVis/full 
M,l 

pVis/full 
M,2 



pVis /full 
M,0 ' 



pvis/full 4a vis/full 



1,1 



pVis/full _8a_ vis/full 



and 



pvis/full pvis/full 



(15) 
(16) 

(17) 
(18) 

(19) 
(20) 



for «• = 0,1,2. 

The cross section predictions avis(£>*^) depend not only on the kinematics of the charm quark 
production mechanism but also on the fragmentation of the charm quark into Z)*^ mesons. The 
charm quark fragmentation function to D*^ mesons has been measured by HI Il32l using the pro- 
duction of D*-mesons with and without associated jets in DIS. In the calculation of avis(£>*^) the 
longitudinal fragmentation is performed in the y* — p centre-of-mass frame, using the fragmen- 
tation function of Kartvelishvili et al. [[33l which is controlled by a single parameter, a^:. This 
parameter has been determined for two different regions of the partonic centre-of-mass energy 
squared, s, depending on the jet requirements made in the different analyses. The fragmentation is 
observed to become softer with increasing s as expected from the perturbative evolution of the frag- 
mentation function not implemented in HVQDIS. This is accounted for in the current analysis by 
using different values of a^: corresponding to the measurements for two ranges of s, as listed in ta- 
bled] The limits on the ranges in s are determined with HVQDIS by applying the jet requirements 
of the individual analysis on parton level. The Uk parameters and the s ranges are varied according 
to their uncertainties to evaluate the corresponding uncertainty on a^-,,. The charmed hadrons also 



s range [GeV-^] 




measurement 


s< (70 ±40) 


6.1 ±0.9 


[32] D*, DIS, no-jet sample 


(70±40) <5 


3.3±0.4 


[l3a£>*,DIS,jet sample 



Table 1 : The parameter aK of the Kartvelishvili et al. fragmentation function for the D* mesons employed 
in our analysis versus parton centre-of-mass energy squared s, equation 

receive a transverse momentum, kj, with respect to the charm quark direction according to 

f{kT)=kT0Xp{-2~kT/{kT))- (21) 

The transverse momentum average kr is chosen as 0.35 ±0. 15 GeV^, in line with the experimental 
results on hadron production in collisions ll34l - [39ll . A fragmentation fraction of charm quarks 
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into D mesons of /(c — t- = 0.2287 ±0.0056, is used as determined by averaging the e+e 
and ep results [|40ll . In order to evaluate the fragmentation model uncertainty on ayis (Z)*±)) the 
parameters kj, CLk, and /(c — )■ D*+) are varied within the uncertainties quoted above and each 
variation is considered as a correlated uncertainty in £vis. In a few cases, a symmetric variation 
of the model parameters results in an asymmetric uncertainty on the cross section. In such cases, 
the largest absolute value of the uncertainty is assigned. The dominant systematic uncertainty is 
arising from the variation of the fragmentation function parameter a^:. 

The measurement p6l used in the fit is based on determination of the vertex displacement of 
the tracks and covers essentially the full phase space. However, the extrapolation to the cross 
sections of charm and beauty quark production is performed with a Monte-Carlo simulation and 
the determination of a correction at NLO as in case of D*^ analysis is not possible. We assume 
this correction to be small and no extrapolation factor is applied to these data. 



4 Results 

The NLO variant of our analysis applies 0{a^) corrections to the heavy quark electro-production 
Wilson coefficients ^ and the NLO ABMl 1 PDFs. For the NNLO variant of our analysis we use 
the NNLO ABMll PDFs and apply the corrections up to 0{a^) as discussed in section [H Since 

the NNLO Wilson coefficients / in equation ([8]) are still approximate and affected by a residual 

uncertainty parameterised by c^-'^ and c^^'^ we define a particular shape for c^^ as a linear 
combination of the two envelopes with a parameter dpj interpolating between these two options, 

4J = {i-d^)4f^+d^4f\ ill) 

The consistency of the prediction using results of the fit with the data of ||4l[5]| on o^^^^^^^ = 
<3^j\^{D*^)/ f{c -^■ D*^) and of (H on o^gjj is illustrated in figures [T] and [2l respectively. In the 
case of NLO, the fit quality is very good for each data set with the total value of = 60 for the 
number of data points (NDP) equal to 60. The consistency of the NNLO prediction using the fit 
results with the data is equally good as for the NLO variant of the fit (cf. figures [T] and |2l), with the 
best fit value of y}/NDP = 63 /60 being achieved for = -0.6. The variant of the NNLO Wilson 

coefficient given by C2 /' (equations (4.18) and (4.22) in [[T6l ) is clearly disfavoured by the data 

7 (2) A 

with % = 156, while the description provided by the variant /' (equations (4.17) and (4.21) 



in [[T6l ) with %^ = 72 is comparable to the best one. Therefore we assign to the value of mdmc) at 
NNLO a conservative theoretical uncertainty of 100 MeV which accounts for the incomplete cur- 
rent knowledge of the NNLO Wilson coefficient and which corresponds to the difference between 
the values of m^.(m^.) obtained with dp^ = —0.6 and dp^ = 0. 

As a result of the QCD analysis, we have determined the MS charm quark mass at NLO and at 
approximate NNLO accuracy in QCD, 

m,.(m,.) = 1.27 ±0.05(exp)+jJ;j!f (scale) NLO, (23) 

m,(m,) = 1.36±0.04(exp)l!;;!!J(scale) ±0.10(th) NNLOapprox- (24) 

The experimental uncertainty is obtained from the propagation of all uncertainties in the data with 
account of their correlations. The fragmentation model uncertainties relevant for the data [|4l|5|| are 
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also included on the same footing. To estimate the influence of the PDF uncertainty on the preci- 
sion of the mc(mc) determination, the NLO analysis is repeated for each of the ABMll PDF set 
members, representing the la uncertainty in the fitted PDF parameter^. The differences between 
the values mc{mc) obtained in this way and the ones obtained with the central PDF member are 
added in quadratures. The resulting uncertainty on mc(mc) is about 10 MeV which is very small 
in comparison to the other uncertainties. 

As the fit is performed at fixed order in perturbation theory, one needs to account for missing 
contributions beyond the order considered. This uncertainty is calculated in the standard manner 
from the variation of the renormalisation and factorisation scales within the range given in equa- 
tion d?]). The scale variation is performed simultaneously in equation (|3]) and in the calculation of 
the extrapolation factors £vis. For the NNLO case the parameter rf^y is fixed at the value of —0.6, 
which is preferred by the fit. Due to the sensitivity of Eyis on the scale variation, the resulting 
uncertainty on mc(mc) does not improve significantly from the NLO prediction to the NNLO one. 
In case of the NNLO result, the theory error on mc(mc) of 100 MeV due to lacking knowledge on 
the NNLO Wilson coefficient is significantly larger than the uncertainty due to scale variation. The 
necessary theory computations to remedy this unsatisfactory situation are discussed in [fT6l . 

The NNLO predictions for charm electro-production obtained for the two variants of the fit are 
shown in figure |3] for the kinematics of the HERA collider experiments. The two variants of the fit 
are not very different at large which illustrates the moderate potential of the data [|4]-[6l to sepa- 
rate between different shapes of the NNLO term in the heavy-quark Wilson coefficients. However, 
at small the predictions diverge therefore the combined HI and ZEUS data may provide an 
additional constraint on this term and can improve the accuracy of the NNLO determination of nic. 

In the previous analysis [|3l[8l, values of the charm quark MS mass of mc{mc) = 1.26 ± 
0.09(exp) ±0.11(th) at NLO and mc(mc) = 1-01 ±0.09(exp) ±0.03(th) at NNLOapprox was ob- 
tained. While the NLO results of the present and the previous analyses are consistent, the present 
result at NNLOapprox is significantly larger. This shift in the central value is due to several sources. 
First, in the present study we use a significantly improved theory of heavy-quark production in 
DIS. As mentioned above, we reconstruct the three-loop Wilson coefficient C2 / based on all avail- 
able information in various kinematic limits lfT6ll . in contrast to earlier studies |l3l[8l which relied 
only on the threshold approximation Cj^j'*^. Also, the analysis of f3lf8]| was performed as a variant 
of the ABKM09 global analysis [[HI of world DIS data, to which the HERA experiments, at that 
time, only contributed the inclusive data of [I411I42II . The present determination of mc{mc) is based 
on the ABMl 1 fit [[3, which incorporates the new combination of the HERA run I inclusive DIS 
data from the HI and ZEUS experiments [[431 and to which we add the new precise data sets on 
DIS charm production [[4]-[6l discussed at above. 

The values of mc{mc) in equations ([23l) and ([24l) determined at NLO and at NNLOapprox are 
close to the world average of mc{mc) = 1.275 ±0.025 GeV as evaluated by PDG [[2l and illustrated 
in figure HI The individual charm-quark mass determination used in the averaging procedure0 
are summarised in table [2l They are based on different theoretical methods but are essentially 
limited to two approaches only. They have either been obtained non-perturbatively with the help 
of lattice QCD simulations with a given number of dynamical fermions. Or, in case perturbative 

^ The parameterisation of the proton structure functions in the ABMll PDF set |7| incorporates the effect from 
higher twist terms described by operators of dimension six in the framework of Wilson's operator product expansion. 

^Note that the PDG converts all results to mc{mc) using scheme transformation ( [T2] | at two loops in QCD together 
with the value a, {mc) = 0.38 ±0.03. 
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mcinic) 


reference 


determination method 


1.261 ±0.016 


NARISON mi 


QCD sum rules for vector current correlator 


1.278 ± 0.009 


BODENSTEIN gU 


QCD sum rules for vector current correlator 


1 90 +0.07 
^■^'^ -0.06 


LASCHKA pel 


lattice QCD; charmonium spectrum 


1.196 ±0.059 ±0.050 


AUBERT [471 


inclusive spectra in semileptonic B-decays 


1.28 ±0.04 


BLOSSIER 148] 


lattice QCD («/ = 2); hadron spectrum 


1.273 ± 0.006 


MCNEILE ||49J 


lattice QCD (nf = 2+1); pseudo-scalar current correlator 


1.279 ±0.013 


CHETYRKIN H 


e+e^ cc cross-section and QCD sum rules. 


1.25 ±0.04 


SIGNER El 


non-relativistic QCD sum rules and e+e^ cc cross- 
section near threshold. 


1.295 ±0.015 


BOUGHEZAL 121 


e+e^ — > cc cross-section. 


1.24 ±0.09 


BUCHMULLER El 


global fit to inclusive B-decay spectra. 


1.224 ±0.017 ±0.054 


HOANG 153 


global fit to inclusive B-decay data. 



Table 2: Individual results of the charm-quark mass determinations used to compute the world average in 
the MS scheme H. 



QCD predictions have been applied, the extractions of mc{mc) listed in table |2] are limited to 
processes with time-like kinematics, i.e., cross-section data from e+e -collisions exposed to QCD 
sum rule analyses or data on 5-decays. The latter method carries also an intrinsic dependence on 
the uncertainty of the Z?-quark mass, while the former, i.e., the QCD sum rule analyses fix the value 
of as{Mz) to the world average including the associated very small erroiE]. It has been shown ll55l . 
that the systematic shift of mc{mc) due to the value of as{Mz) in QCD sum rule analyses is quite 
sizable. 

5 Conclusions 

The high precision measurements of charm quark production in DIS at HERA offer an attractive 
possibility to extract the charm-quark mass in the theoretically well-founded MS scheme. The 
resulting experimental precision of the present determination based on DIS data is substantially 
improved with respect to previous analyses. It is now compatible with the theoretical uncertainty 
when comparing to the QCD predictions at NLO and it is significantly smaller than in case of the 
approximate NNLO QCD predictions. The latter suffer from missing information on the three- 
loop Wilson coefficients at small-jc and small values of and imply an additional theoretical 
uncertainty on /n^ of 100 MeV. 

In comparison to the measurements used to define the world average up to now, our results 
add complementary information from scattering processes with space-like kinematics and provide 
an important test of the QCD dynamics. The kinematic range covered by the DIS data allows for 
the extraction of lUc in the MS scheme well within the regime of validity of perturbative QCD. 
This analysis accounts for the full correlation of the dependence on iric with other non-perturbative 

^ The current world average of (Xs{Mz) is the result of an arithmetic average of high precision determinations at 
least to NNLO which are only marginally compatible within their quoted errors, see e.g., 12117]. 
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parameters, most prominently the value of strong coupling constant as{Mz) and the gluon PDF. 
Future improvements on the accuracy of the mc{mc) extractions from DIS data rely mostly on the 
theoretical progress for the three-loop Wilson coefficients especially at small x. Also the experi- 
mental uncertainty can still be reduced by combining the charm production data from the HI and 
the ZEUS collaborations. 

In summary, the charm-quark mass extractions from DIS data may challenge the accuracy of 
QCD sum rules analyses once the uncertainties on all non-perturbative parameters are accounted 
for on equal footing. 



Acknowledgements 

This work has been supported in part by Helmholtz Gemeinschaft under contract VH-HA-101 (Al- 
liance Physics at the Terascale), VH-NG-401 {Young Investigator group "Physics of gluons and 
heavy quarks "), by the Deutsche Forschungsgemeinschaft in Sonderforschungsbereich/Transregio 9 
and by the European Commission through contract PITN-GA-2010-264564 (LHCPhenoNet). 



References 

1. 1. Bigi, M. A. Shifman, N. Uraltsev, and A. Vainshtein, Phys.Rev. D50, 2234 (1994), [hep^h/940236q 
Particle Data Group, J. Beringer ef aL, Phys. Rev. D86, 010001 (2012). 
S. Alekhin and S. Moch, Phys.Lett. B699, 345 (2011). larXiv:1011.5790[ 
HI Collaboration, F. Aaron et ai, Phys.Lett. B686, 91 (2009). larXiv:0911. 39891 
HI Collaboration, F. Aaron ef a/., Eur.Phys.J. C71, 1769 (201 D. larXiv: 1106.10281 
HI Collaboration, F Aaron et ai, EurPhys.J. C65, 89 (2010) jarXiv:0907.2643l 
S. Alekhin, J. Blumlein, and S. Moch, (2012'). larXiv:1202.2281l 
S. Alekhin and S. Moch, (201 D. larXiv: 1107.04691 

E. Laenen, S. Riemersma, J. Smith, and W. van Neerven, Nucl.Phys. B392, 162 (1993). 



S. Riemersma, J. Smith, and W. van Neerven, Phys.Lett. B347, 143 (1995), hep-ph/94 11431 



M. Buza et ai, Nucl.Phys. B472, 611 (1996), hep-ph/9601302 



I. Bierenbaum, J. Blumlein, and S. Klein, Nucl.Phys. B780, 40 (2007), [hep^h/0703285| 

M. Buza, Y. Matiounine, J. Smith, and W. van Neerven, Eur.Phys.J. CI, 301 (1998), [h^h/9612398| 



S. Alekhin, J. Blumlein, S. Klein, and S. Moch, Phys.Rev. D81, 014032 (2010). larXiv:0908.27661 
S. Forte, E. Laenen, P Nason, and J. Rojo, Nucl.Phys. B834, 116 (2010). larXiv: 1001.23 121 
H. Kawamura, N. L. Presti, S. Moch, and A. Vogt, (20 1 2). larXiv: 1205.57271 
N. Gray, D. J. Broadhurst, W. Grafe, and K. Schilcher, Z.Phys. C48, 673 (1990). 



K. Chetyrkin and M. Steinhauser, Nucl.Phys. B573, 617 (2000), [hep^h/991 1434 

K. Melnikov and T. v. Ritbergen, Phys.Lett. B482, 99 (2000),'hep-p h/9912391| 

U. Langenfeld, S. Moch, and P Uwer, Phys.Rev. D80, 054009 (2009). larXiv:0906.5273[ 

M. AUev et ai, Comput.Phys.Commun. 182, 1034 (201 1), arXiv: 1007.13271 



J. Vermaseren, S. Larin, and T. van Ritbergen, Phys.Lett. B405, 327 (1997), hep-ph/9703284 



10 



K. Chetyrkin, Phys.Lett. B404, 161 (1997), [hep^h/9703278 



B. Harris and J. Smith, Nucl.Phys. B452, 109 (1995), |hep^h/9503484| 
E. Laenen and S. Moch, Phys.Rev. D59, 034027 (1999), |hep^h/9809550| 



S. Alekhin and S. Moch, Phys.Lett. B672, 166 (2009). larXiv:0811.1412l 

I. Bierenbaum, J. Bliimlein, S. Klein, and C. Schneider, Nucl.Phys. B803, 1 (2008). EXiv:0803.0273l 

I. Bierenbaum, J. Bliimlein, and S. Klein, Nucl.Phys. B820, 417 (2009'). larXiv:0904.3563l 

I. Bierenbaum, J. Bliimlein, and S. Klein, Phys.Lett. B672, 401 (2009'), larXiv:0901.0669l 

J. AbUnger et al, Nucl.Phys. B844, 26 (201 n. larXiv: 1008.33471 

S. Catani, M. Ciafaloni, and F. Hautmann, Nucl.Phys. B366, 135 (1991). 

HI Collaboration, F. Aaron et al, EurPhys.J. C59, 589 (2009). larXiv:0808.1003l 

V. Kartvelishvih, A. Likhoded, and V. Petrov, Phys.Lett. B78, 615 (1978). 

ALEPH Collaboration, R. Barate et al, Phys. Rep. 294, 1 (1998). 

DELPHI Collaboration, P Abreu et al, Z.Phys. C73, 1 1 (1996). 

MARCH Collaboration, G. S. Abrams et al, Phys. Rev. Lett. 64, 1334 (1990). 

PLUTO Collaboration, C. Berger et al, Phys. Lett. B82, 449 (1979). 

PLUTO Collaboration, C. Berger et al, Phys. Lett. B97, 459 (1980). 

TASSO Collaboration, R. Brandelik et al, Phys. Lett. B94 (1980). 

E. Lohrmann, (201 D. larXiv: 11 12.37571 



HI Collaboration, C. Adloff ef aL, EurPhys.J. C21, 33 (2001), [hep^x/0012053 

ZEUS Collaboration, S. Chekanov et al, EurPhys.J. C21, 443 (2001), [hep^x/0105090| 

HI and ZEUS CoUaboration, F Aaron et al, JHEP 1001, 109 (20 10). EXiv:09 11.08841 

S. Narison, Phys.Lett. B706, 412 (2012). |arXiv: 1105.29221 

S. Bodenstein et al, Phys.Rev. D83, 074014 (201 D. larXiv: 1102.38351 

A. Laschka, N. Kaiser, and W. Weise, Phys.Rev. D83, 094002 (201 D. EXiv: 1102.09451 

BABAR Collaboration, B. Aubert et al, Phys.Rev. D81, 032003 (2010). larXiv:0908.0415l 

ETM Collaboration, B. Blossier et al, Phys.Rev. D82, 1 14513 (2010), larXiv: 1010.36591 

C. McNeile etal, Phys.Rev. D82, 034512 (2010). EXiv: 1004.42851 

K. Chetyrkin et al, Phys.Rev. D80, 074010 (2009). EXiv:0907.2110l 

A. Signer, Phys.Lett. B672, 333 (2009). larXiv:0810.1152l 

R. Boughezal, M. Czakon, and T. Schutzmeier, Phys.Rev. D74, 074006 (2006), |hep^h/0605023[ 
O. Buchmuller and H. Fliicher, Phys.Rev. D73, 073008 (2006), [hep^h/0507253| 



A. H. Hoang and A. V. Manohar, Phys.Lett. B633, 526 (2006), [hep^h/0509195[ 



B. Dehnadi, A. H. Hoang, V. Mateu, and S. Zebarjad, (201 D. larXiv: 1102.22641 



11 



cc,VIS 




Figure 1: The visible reduced charm quark production cross section, as determined from the 
measurement of D* meson production by the HI collaboration in comparison to the predicted 
cross sections using the NLO (solid line) and NNLO (dashed line) variants of this analysis. 
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Figure 2: The reduced charm quark production cross section as measured by the HI collaboration 
using the vertex information of inclusive track production in comparison to the predicted cross 
sections using the NLO (solid line) and NNLO (dashed line) variants of this analysis. 
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Figure 3: The NNLO predictions for the reduced cross section of the charm electro-production as 
a function of x for different Q^, for kinematics of the HERA collider experiments. The solid line 
represents the variant of the fit using = —0.6 and mc{mc) = 1 .36 GeV. The result corresponding 
to the fit using dp^ = and mc{mc) = 1.27 GeV is shown by a dashed line. 
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Figure 4: The charm quark MS mass as determined at NLO (open square) and at NNLOapprox 
(filled square) using DIS data including the measurements of open charm production at HI. The 
inner (outer) error bars represent experimental (total) uncertainties. To obtain the total uncertainty, 
the experimental and the theory uncertainties arising due to variation of the scales are added in 
quadrature. For comparison, the world average evaluated by the PDG f2| is shown by the shaded 
band. The measurements entering the world average determination are represented by open trian- 
gles. 
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